Luttinger Liquid State with Effective Attractive Hard-Core Interaction 
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An exact solvable 'zig-zag' ladder model of degenerated spinless fermions is proposed and solved 
exactly by the means of the Bethe ansatz. An effective attractive hard-core interaction and direct 
Coulomb repulsion of fermions on the nearest-neighbor sites of different chains induce a new phase 
state of the ladder. We give a detailed analysis of the exact phase diagram at zero temperature, 
that is characterized by two phases at the filling exceeding 2/3: itinerant and 'frozen' fermions, for 
which two different species reside in spatially separate regions. The critical exponents describing 
asymptotic behavior of the correlation functions are calculated using the Bethe ansatz and conformal 
field theory. It is shown also, that the density of the magnetization of the corresponding zig-zag 
spin- 1/2 ladder has a jump equal to 1/3 at the magnetic field equal to zero. 
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The metallic single-walled carbon nanotubes 
(SWNT's) represent one-dimensional (ID) systems 
(in the sense of their electronic properties), in which 
one may expect to observe physical phenomena char- 
acteristic for strong electron correlations. Informer 
studies attention has been mostly concentrated on the 
realization of the Tomonaga-Luttinger liquid state in 
these systems. This conclusion follows from the analysis 
of their transport propertiesi*^ and from the direct 
observation of electronic states near the Fermi energy 
in the high-resolution photoemission experiments^. The 
large value of the critical exponent 6 (6 = 0.3 — O.Sii^i^), 
extracted from the experimental results, which describes 
the asymptotic behavior of the spectral function near 
the Fermi energy, indicates that large density-density 
correlations dominate in the system. In other words, 
SWNT's are ID electronic systems with strong repulsive 
interaction between fermions. The strong correlations 
work wonders in such electron liquid, therefore a tra- 
ditional Luttinger liquid approach^ that is relevant in 
the case of a weak electron-electron interaction can not 
explain the mechanism of interactions inherent in carbon 
nanotubes^. It is necessary to use exactly solvable 
models or numerical calculations of many-body fermion 
systems with strong interactions for the complete un- 
derstanding of such intriguing behavior of SWNT's. We 
shall exploit this conception and develop an approach 
based on the exact solvable zig-zag ladder model of 
spinless fermions. An effective attractive hard-core and 
direct repulsive interactions are intrinsic in the model 
considered. As it will be shown below, the interplay 
of these interactions gives a rise to the phase diagram, 
which includes interesting strongly correlated states of 
fermions at the high filling. The values of Q calculated in 
the framework of the model proposed are in qualitative 
agreement with Refiii^i^. 

The metallic state of SWNT's is realized in the two 
possible high symmetry structures for nanotubes, known 
as zig-zag and 'armchair' (m = for all zig-zag tubes 
while n = m for all armchair tubes^). The coupled 
chains in the form of the ladders are hard lattice mod- 



els of strongly correlated systems exactly solvable in 
uncommon cases. Indeed, such systems are frequently 
frustrated^. The models of the coupled spin chains® and 
spin ladders have been studied intensively^ : ' ^ ^ . Never- 
theless, there are still many open questions on the gaps 
in the spectrum of collective excitations and the behavior 
of these systems away from half-filling. In this Letter we 
consider the model for zig-zag nanotubes in the form of 
the zig-zag ladder of degenerated spinless fermions taking 
into account the hopping of fermions along the chains and 
the interaction between fermions on the nearest-neighbor 
lattice sites of different chains. The model Hamiltonian 
is given by 

eiMi + e2M2 + ^Y^^rij + n,+i)n^.+ i , (1) 
j 

where Cj and are operators of spinless fermions at site 
j shifted on a half-lattice constant for different chains, ac- 
cordingly the first and the second chains, the hopping in- 
tegral is equal to unity, ei and £2 define a relative shift of 
the fermion subbands of the chains, J is the repulsion in- 
teraction of fermions on the nearest-neighbor sites of dif- 
ferent chains, the particle number operator for fermions 
is defined by nj = CjCj. The summation extends over all 
sites of the chains of the length L, we assume periodic 
closure. The hopping of fermions through the sites oc- 
cupied by fermions of another chain is forbidden in the 
hopping term of the Hamiltonian. The Hamiltonian con- 
serves the total number of fermions of each chain Mi and 
M2, N = Ml -I- M2 is the total number of particles. At 
ei — e2 the chains are identical and the Hamiltonian (1) 
is reduced to the zig-zag spin-1/2 ladder in the magnetic 
field if - ei n - -EJ(^J^J+l+'5j'5jVl)(5|+i + l/2)- 
^ Ej + iSj+S^+i)S^^,. By way of illustration 
of the structure of the zig-zag ladder and interaction the 
chains are shown in Fig. 1. 

We now turn to the diagonalization of the model 
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FIG. 1: Graphical representation of the zig-zag ladder; and 
1 denote the value of the hopping integral along the chains 
between the nearest-neighbor sites at given configuration of 
fermions. 



limit the rapidities Xj are closed spaced and may be re- 
garded as continuous variable. The distribution functions 
of the charge rapidities pi,2(A) for each chain are defined 
through the integral equations of the Fredholm type 



Pl,2(A) 



dA'i?2(A- A')p2a(A') = 
(l + m2,i/2)i?i(A), 



(4) 



Hamiltonian by the coordinate Bethe ansatz formalism. 
The Schrodinger equation is solved with the ansatz 



'tp{xi,X2, ...,xn) = y^(-l)^^(P)exp 



N 



,(2) 



where the P summation extends over all the permuta- 
tions of the momenta {kj} of the particles. 

According to the form of the interaction in (1) the 
fermions of one chain are scattered on the fermions of 
another only, the two-particle scattering matrix of the 
spinless fermions with momenta hi and kj is given by 

l-l-exp[8(fci+fcj)] + J exp(zfcj) 



Sij 



exp[2(A:j kj)] i+exp[i(ki+kj)]+Jexpiiki) ■ 



ditional constrain on the form of the interaction in (1) 
follows from the calculations of the configurations with 
three particles in different chains on the neighboring 
lattice sites (it is sketched in Fig. 1). Direct calcula- 
tions show that if the model Hamiltonian (1) does not 
contain a direct interaction of fermions on the nearest- 
neighbor sites of the same chain (when the site between 
these particles is occupied), the relation between ampli- 
tudes of the three-particle wave function holds for arbi- 
trary configuration of the particles A(kp'i, kp>2, kp's) = 
S{kpj,kpj+i)A{kpi,kp2,kp3) (here P is an arbitrary 
permutation and P' = P{j,j + 1)). 

The energy eigenstates are characterized by sets of the 
charge rapidities = 1, ...,N) for the particles, that 

satisfy the following Bethe equations written below in the 
case of a weak interaction J/2 = cos rj < 1 for the charge 
rapidities of the 1-st chain {{Xj}, j = 1, Mi) 



sinh i {Xj 



■irj) 



L-t-iM2 



sinh ^{Xj — if]) 

^ sinh i(Aj - Aj -f 2ir]) 
l\smh^{Xj-Xi-2ir,)- 



(3) 



where exp(ifcj) 



sinh ^(Aj-t-iry) 



> '^2 = ki is the mo- 



sinh —iri) ' 

mentum of the 2-nd chain. The eigenenergy of the ladder 
is defined hy E = -2X;jLi coskj + eiMi -|- eaMa. The 
charge rapidities of the 2-nd chain are solutions of the 
similar Bethe equations. 

The fillings of the fermionic subbands are defined by 
the chemical potential of the system, which is the same 
for different subbands at their partial filling, thus Mi = 
M2 at ei = 62 for identical chains. In the thermodynamic 



1 sin(n?7) 
27r cosh A— cos(n77) ' 



= f^Api,2(A) are the 



with the kernel RnW = 

TO12 = Mi^2/L 
densities of fermions in the chains, £1.2 = 
nil, 2^1, 2 — 47rsin77 J^)^^ ^ dXRi{X)pi^2{^) are the densities 
of the ground-state energy of the chains. The similar 
Bethe equations take place in the case of a strong 
repulsive interaction at J/2 = cosh/x > 1. The corre- 
sponding integral equations are defined by the kernel 

P (\'\ — 1 sinh(Ttjj) 
-n-nl^j — 271- cosh(n/i)-cosA- 

The low-density region of the phase diagram, when the 
distance between particles is large compared to the hard- 
core radius, describes itinerant fermions in the Coulomb 
potential J, since an effective force of the the hard-core 
interaction is insignificant. Complicated physical picture 
arises in region of high densities of fermions, since radi- 
cally different situation is realized in the model at high 
filling of the chains when the average distance between 
particles is compared with the hard-core radius and the 
hard-core interaction dominates. The integral equations 
(4) are defined for the density of particles less than a 
'half-filling' equaled to 2/3 for identical chains. For in- 
vestigation of fermionic states of the ladder at larger 
filling let us use the combined electron-hole symmetry 
Cj (—!)■' c] for the Hamiltonian (1) and obtain H ^ H' 

n' = J2[-c]cj+i - c]_^_T^Cj + ^J{nj + nj+i)]nj+i 
j 

-JN + ei{L-Mi) + e2{L-M2). (5) 

The Hamiltonian (5) describes the hole states in the lad- 
der, whereas the fermion states of the first and hole states 
of the second chains are defined by the following Hamil- 
tonian 



n' = - ^(c]c,-+i + c]+ic,)n^-+ 1 + (ei + J)Mi - 



J2 ic]cj+i + cj+iCj)(l - rij^i) + 



j=l/2 



(£2 - J){L - M2) - ^JY^iuj + nj+i)nj+i. (6) 



The phase state of the model is result of the interplay 

of the repulsive Coulomb interaction J and effective at- 
tractive hard-core interaction with a negative hard-core 
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FIG. 2: Ground-state energy as a function of the density of 
spinless fermions. Individual curves are labelled by value of 
J=0;0.5;l;2;3;5;10, a dotted curve corresponds to free fermion 
state. 



partial densities of fermions 



1.0 
0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 



I I I I 1 I I I I 1 I I I I 1 I 

one component state 


mixed state/ : 


Itinerant fermions / 




/ frozen fermions 

/ \ 





0.0 0.2 0.4 0.6 0.8 1.0 
density of fermions 

FIG. 3: Phase diagram in the coordinates of the densities of 
itinerant and frozen fermions as a function of the density of 
fermions. 



radius equaled to -1/4. We shall show, that phase sep- 
aration occurs at high densities of particles [n > 2/3). 
First of all let us answer on the question: what is it 
the effective attractive hard-core interaction? A repul- 
sive hard-core interaction forbids two particles being at 
distances less or equal to doubled hard-core radius I, that 
corresponds to equaled to zero two-particle wave function 
ilj{xi,Xj) at \xj — Xi\ < I. The interaction in the Hamil- 
tonian (1) forbids the tunneling of particles of different 
chains on distance of a haft-lattice constant, other words 
the two-particle wave function is equal to zero after the 
result of the hopping (or the current is equal to zero), 
then ip{x,x) ^ automatically in the case for a nontriv- 
ial solution for the wave function. Thus this condition 
on the two-particle wave function leads to an effective 
attraction between two particles of different chains with 



the same coordinates in contrast to a traditional repul- 
sive hard-core interaction for that ip{xj,Xj) — 0. 

Using exact solutions of the Hamiltonians (5), (6) we 
can calculate the phase diagram at arbitrary fillings of 
the chains, below we consider identical chains. The 
Hamiltonians (5) is defined the phase state of fermions 
with frozen hopping. The Bethe function (2) is the 
solution of the Schrodinger equations with equaled to 
zero energy and the two-particle scattering matrix S'^j ~ 
exp[|(fcj — ki)]. Really this phase corresponds to full- 
filled states of the chains, other words these solutions 
describe the phase separation of the chains at n > 2/3: 
the full-filled states of fermions in the chains with frozen 
hoppings and the 'old' phase of itinerant fermions defined 
by the Bethe equations (3) on the restricted chains. The 
density of the ground-state energy is plotted in Fig. 2 
at ei = 62 = 0. At n > 2/3 the density of 'frozen' 
fermions is equal to Sn = 3n — 2, the corresponding den- 
sity of the itinerant fermions An = n — 6n — 2(1 — n) 
(see Fig. 3). In a separate point J — > oo the scattering 
matrix is defined by conditions '!/'(a^j±i , Xj) = 0, that cor- 
responds to a constrain on the tree-particle wave function 
'0(a;j+i, a;^^! , Xj) — also, therefore the maximal filling 
of chains is equal to 1/3. In the case of a strong repulsive 
interaction J > 2 both phases are dielectric, the phase 
of itinerant fermions is metallic for a weak interaction 
J < 2. 

The Luttinger liquid parameter K^J^ is defined by the 



X?,2(A). 



The inte- 



dressed charge Xi,2(A) via K^p'"^ 
gral equations for the dressed charge functions Xi,2(^) is 
distinguished from equations (4) by the driving term 



Xi,2(A) 



dA'i?2(A - A')x2,i(A') - (1 + m2,i/2) .(7) 
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FIG. 4: The critical exponents Kp and at J=0;0.5,l;2;3;10. 



According to (7) the repulsive interaction J decreases 
the value of Kp. In contrast to traditional hard-core 
interaction, that decreases the value of Kp extremely^, 
an effective attractive hard-core interaction increases Kp. 

-> corresponds to free spin- 
In the case of the identical 



The low density limit mi 2 
less fermions with Kn 1 
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FIG. 5: Density of the magnetization of the spin ladder model 
as a function of the magnetic field at J=0.01; 1; 2; 2.5; 3. 

chains the system of integral equations (7) is reduced to 
one and we obtain the solution for Kg = |— ^ at half 
filling mi 2 = 2/3 using the Wiener-Hopf method. The 
maximal value of Kp = 16/9 is realized in the J ^ 
limit, that corresponds to © = 0.17. The critical expo- 
nent Q is related to the Luttingcr parameter Kp through 
8 = {Kp — 2 + l/Kp)/2. Whereas, the minimal value 
Kp = 4/9 is obtained in the J — > oo limit, we have twice 
larger vahic of the critical exponent = 0.35 in this 
case. Numerically solving integral equations (4), (7) we 
calculate the critical exponents Kp and Q for each chain 
as function of the density of spinless fermions for differ- 
ent values of J. Numerical calculations of the exponents 
are shown in Fig. 4 for illustration of the behavior of the 
critical exponents for arbitrary densities of particles and 
interaction. In an effective attractive regime at a weak 
interaction (J < 2) the value of Kp decreases from 16/9 
to 1 with increasing J from to J = Jq (the value of Kp 
depends of the density fermions therefore Jq depends of 
n also, according to simulations 1.7 < Jq < 1.9), whereas 



the value of O decreases from 0.17 to zero. The phase of 
itinerant fermions is metallic at n > 2/3, therefore the 
critical exponents independent of the density fermions at 
n > 2/3. For J > Jq Kp decreases with n and reaches 
a minimal value 4/9 at half-filling, the value of Kp de- 
creases with J also. The value of 8 increases and reaches 
a maximum value 0.37 at half-filling, both phases are di- 
electric for J > 2 at n > 2/3. At half-filling the value 
J = 2 separates the phase of itinerant fermions on gapless 
at J < 2 and gaped at J > 2. The behavior of the den- 
sity of the magnetization as a function of the magnetic 
field, calculated for difi^erent values of the interaction J, 
is shown in Fig. 5. The magnetization of the system has 
a jump at H=0 equal to 1/3. This jump corresponds to 
the phase of frozen fermions with the chemical potential 
equal to zero. 

In summary, we presented the exact phase diagram of 
the zig-zag ladder model for arbitrary Coulomb repulsion 
interaction and density of fermions obtained on the basis 
of the exact solution of the model. It is shown that intrin- 
sic effective attractive hard-core interaction and Coulomb 
repulsion form the Luttinger liquid state in the ladder. 
Due to special role of the effective attractive hard-core in- 
teraction, two phases of the itinerant and frozen fermions 
are separated at the high density, when filling exceeds 
2/3. As a result, the system exhibits two types of behav- 
ior in high density region: either band insulator in the 
case of the strong interaction at J > 2 or the Luttinger 
liquid state for one phase and band insulator for another 
at J < 2. The values of bulk exponents calculated in 
the framework of the model adopted are in a reasonable 
agreement with the measurements in SWNT's. 
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